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Introduction 



Let us consider the propagation of scalar waves with the speed v > given by the 
^ . wave equation Uu — v'^Au = outside a compact domain D in the Euchdean space 

=3 : K'^. Let us put Q = R'^\D. We can assume for example Neumann boundary 

conditions. We will denote by Af^ the (self-adjoint) Laplace operator with the 
boundary conditions. So our stationary wave equation is the Helmholtz equation 

v'Anf + u;'f = (1) 

^ ■ with the boundary conditions. We consider a bounded interval / = C 

rn 



]0, +oo[ and the Hilbert subspace T-Lj of L (fi) which is the image of the spectral 



I projector Pj of our operator —v A 



Let us compute the integral kernel n/(a;, y) of Pj defined by: 

Pif{x)= [ Uj{x,y)fiy)\dy\ 
Jn 

into 2 different ways: 

1. From general spectral theory 

2. From scattering theory. 

Taking the derivatives of Ilj{x,y) w.r. to cj+, we get a simple general and ex- 
act relation between the correlation of scattered waves and the Green's function 
confirming the calculations from [2] in the case where D is a disk. See also [1] 
for other relations between correlations and Green's functions used in passive 
imaging. 
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1 n/(x, y) from spectral theory 

Using the resolvent kernel (Green's function) G{u,x,y) = [{u^ + v'^ An)^^]{x , y) 
for Imu > and the Stone formula, we have: 

n/(a;, y) = Im ( / G{uj + iO, x, y)udu 



TT 

Taking the derivative w.r. to cj+ of n[^^2 ^uj^ji^^ 2/)? we get 

:7-n[^_,<^2](a;, y) = -—lm{G{u + iO, x, y)) . (2) 

2 Short review of scattering theory 

They are many references for scattering theory: for example [3111]. 
Let us define for k G M"^ the plane wave 

eo(x,k) = e'<'^l"> . 

We are looking for solutions 

e{x, k) = eo(a;, k) + e^{x, k) 

of the Helmholtz equation ([T]) in where c^, the scattered wave satisfies the 
so-called Sommerfeld radiation conditionll]: 



e\x, k) = or;™ e~( — , k) + 0(— ) , a; ^ oo . 



\x\ 



The complex function e°°(x,k) is usually called the scattering amplitude. 

It is known that the previous problem admits an unique solution. In more 
physical terms, e(x, k) is the wave generated by the full scattering process from 
the plane wave eo(x,k). Moreover we have a generalized Fourier transform: 



/(x) = (27r)-'^ / /(k)e(x,k)|rfk| 



with 



/(k)= / eiy,k)fiy)\dy\ . 



From the previous generalized Fourier transform, we can get the kernel of any 
function $(— f^A^) as follows: 

m-v'An)]ix,y) = {2n)-' [ <l>{v'k')e{x,k)'^i^\dk\ . (3) 



^As often, we denote k := |k| and k := k/fc 
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3 n/(x, y) from scattering theory 

Using Equation ([3]) with $ = 1/ the characteristic functions of some bounded 
interval / = [w^jW^], we get: 

n,(x, y) = {27^)-" I e(x, k)^(^Mk| . 

J ui- <vk<u) 

Using polar coordinates and defining \da\ as the usual measure on the unit {d — 
1)— dimensional sphere, we get: 

ni{x,y) = (271)-'^ I k'^-^dk I e(x,k)e(i/,k)|rf(T| . 

J ui—<vk<u) J]s?=k^ 

We will denote by ad-i the total volume of the unit sphere in W^: ctq = 2, cxi = 
27r, 0-2 = 47r, ■ ■ ■ . 

Taking the same derivative as before, we get: 

d uj'^~^ r 

:i-n[^^,^2](x,?/) = (27r)"'^^- / e{x,k)e{y,k)\da\ . 

Let us look at e{x, k) as a random wave with k = u/v fixed. The point-point 
correlation of such a random wave C^^^^{x,y) is given by: 

cr"(^,2/) = -^ 

o 

Then we have: 



CT\^,y) = / eix,k)e{y,k)\da\ 

O'd-l Jvk=ui 



du 



n[.. y) = (27r)-^ ^CT\x, y) . (4) 



4 Correlation of scattered plane waves and Green's 
function: the scalar case 

From Equations (|2]) and (jl]), we get: 



(27r)-'^ ^'"'^'^-^ Cr"(^,l/) = -^Im(G(a; + ^0,x,y)) . 



■yscatt / 

•^r.y; = - 

TT 

Hence, we have 

Theorem 1 For the scalar wave equation Uu — v'^Au = outside a bounded 
domain in M.'^, we have the following expression of the correlation of scattered 
wave of frequency u in terms of the Green's function: 

CT'\x,y) = - ^'^'""^'I' lmiGiu + tO^x^y)) . 

For later use, we put 

r)d+l„d-l 

(^d-i 
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5 The case of elastic waves 



We will consider the elastic wave equation in the domain f2: 

Hu — w^u = 0, 

with self-adjoint boundary conditions. We will assume that, at large distances, 
we have 

Hu = —a Au — b grad divu . 
where a and b are constants: 

/X A + /i 
a = — , = 

P P 

with A, fi the Lame's coefficients and p the density of the medium. We will 
denote vp := \/a + b (resp. Vs '■= \/a) the speeds of the P— (resp. 5— )waves 
near infinity. 

5.1 The case n = R^ 

We want to derive the spectral decomposition of H from the Fourier inversion 
formula. Let us choose, for k 7^ 0, by k, ki, ■ ■ ■ , k^^.i an orthonormal basis of 
M.^ with k = I such that these vectors depends in a measurable way of k. Let 

us introduce Pp = kk* the orthogonal projector onto k and = Yl'jZi kjk* so 
that Pp + Ps = Id. Those projectors correspond respectively to the polarizations 
of P— and S*— waves. 
We have 

Uj{x,y) = {2nr'^J^,^^u'^-'du [vp' J^^.^^e^'^^^-y^ P^da+ 
using the plane waves 

e^(a;,k) = e^'^'k 

and 

we get the formula^: 

n,(x,y) = {27,)-^ j^.^^u'^-'duj (v-/i^,^^ |e^(a:,k))(e«(y,k)|da+ 
Vs' Lsk=. k))(eg,,.(y, k)\da) . 

^We use the "bra-ket" notation of quantum mechanics: |e)(/| is the operator x — > {f\x)e 
where the brackets are hnear w.r. to the second entry and anti-hnear w.r. to the first one 
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5.2 Scattered plane waves 

There exists scattered plane waves 

ep(a;, k) = ep(x, k) + ep(a;, k) 

esjix, k) = esjix, k) + Cgjix, k) 

satisfying the Sommerfeld condition and from which we can deduce the spectral 
decomposition of H. 

5.3 Correlations of scattered plane waves and Green's 
function 

Following the same path as for scalar waves, we get an identity which holds now 
for the full Green's tensor ImG(w + iO, x, y): 

Theorem 2 For the elastic wave equation, we have the following expression of 
the imaginary part of the Green's function in terms of the correlation of the 
scattered S- and P- waves: 

lmG{u + iO,x,y) = -7^'^^"' X^,^^ |ep(x, k))(ep(i/, k)|rfa + 

Ejt; 4fc=^ |e5j(a;, k))(e5j(i/, k)^^) , 

with defined by Equation 

This formula expresses the fact that the correlation of scattered plane waves 
randomized with the appropriate weights {vp'^ versus v^"^) is proportional to the 
Green's tensor. Let us insist on the fact that this true everywhere in Q even in 
the domain where a and b are not constants. 
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